Coexistence of double alternating antiferromagnetic chains 
in (VO)2P207 : NMR study 
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Nuclear magnetic resonance (NMR) of "^^P and '''^V nuclei has been measured in a spin-1/2 
alternating-chain compound (VO)2P207. By analyzing the temperature variation of the "^^P NMR 
spectra, we have found that (VO)2P207 has two independent spin components with different spin- 
gap energies. The spin gaps are determined from the temperature dependence of the shifts at '^^P 
and ^^V sites to be 35 K and 68 K, which are in excellent agreement with those observed in the 
recent inelastic neutron scattering experiments [A.W. Garrett et at, Phys. Rev. Lett. 79, 745 
(1997)]. This suggests that (VO)2P207 is composed of two magnetic subsystems showing distinct 
magnetic excitations, which are associated with the two crystallographically-inequivalent V chains 
running along the b axis. The difference of the spin-gap energies between the chains is attributed to 
the small differences in the V-V distances, which may result in the different exchange alternation 
in each magnetic chain. The exchange interactions in each alternating chain are estimated and 
are discussed based on the empirical relation between the exchange interaction and the interatomic 
distance. 



PACS numbers: 76.60.-k, 75.10.Jm, 75.30.Et 
I. INTRODUCTION 

Low-dimensional quantum-spin systems with antifer- 
romagnetic interactions exhibit various interesting phe- 
nomena of collective nature. Vanadyl pyrophosphate 
(VO)2P207 is a gapped antiferromagnet with spin 1/2 
and had long been known as a good model compound 
of a two-leg antiferromagnetic Heisenberg spin ladder .ETu 
Recently, dispersion of the magnetic excitations was mea- 
sured using the inelastic neutron-scattering technique 
and iS|-|found to be inconsistent with a simple ladder 
model.Q Although an alternating Heisenberg antiferro- 
magnetic chain is proposed as an alternative magnetic 
model of (VO)2P207, the suggested direction of the chain 
is not readily known from the crystal structure and the 
excitation spectrum is not fully compatible with a sim- 
ple alternating spin chain. It is still a controversy what 
magnetic model describes best the magnetic behavior of 

(VO)2P207. 

The ladder interpretation of (VO)2P207 is based on 
the intuitive consideration of the crystal structure which 
is shown schematically in Fig. 0. (VO)2P207 belongs 
to the monoclinic space group P2i and has unit-cell di- 
mensions of a = 7.28 A, 6 = 16,59 A, c = 9.58 A and 
(3 = 89.98° at room temperature.^ The vanadium atoms 
are tetravalent to have spin-1/2 degrees of freedom, each 
of which is coordinated pyramidally by oxygen atoms. 
The VO5 pyramids form pairs in opposite phases in the 
h direction by sharing their edges, and are stacked in 
pairs along the a axis to make up a two-leg ladder. The 
ladders are linked by PO4 tetrahedra and form a three- 



dimensional lattice structure. The proposed direction of 
the alternating chain is along the b axis which is perpen- 
dicular to the ladder axis and the chain includes PO, 
tetrahedra as one of the important exchange pathways 
There are eight inequivalent V and P sites in the unit 
cell, so that the ladders and chains are not all crystallo- 
graphically equivalent. 

The magnetic susceptibility of (VO)2P207 has a broad 
maximum around 80 K and decreases rapidly at low 
temperatures. Johnston et al. have pointed out that 
(VO)2P207 has a spin-singlet ground state owing to the 
ladder configuration of magnetic atoms, although the 
spin-1/2 alternating Heisenberg antiferromagnetic chains 
described by the Hamiltonian 
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accurately reproduces the susceptibility of (VO)2P207.El 
Numerical calculations of the susceptibility of both the 
two-leg ladder and the alternating chain showed that 
both models fit the susceptibility of (VO)2P207 quite 
well, so that the susceptibility measurement alone cannot 
distinguish which model is appropriate.0 The calculations 
also showed that the magnitude of the spin gap is differ- 
ent depending on the models, which can be tested by in- 
elastic neutron scattering (INS) experiments. The subse- 
quent INS measurement on powder sample of (VO)2P207 
has detected a spin gap of 3.7 meV, which favors the lad- 
der model E 

More recently, however, Garrett et al. have presented 
a surprising INS result on the magnetic excitations from 



an array of a lot of small single crystals.a They found 
that the spin excitation is predominantly one dimen- 
sional but is directed along the crystalline b axis rather 
than the ladder (the crystalline a) axis. The exchange 
interaction along the a axis is ferromagnetic and is or- 
der of magnitude smaller than that along the b axis, 
and there is essentially no dispersion along the c axis. 
As the most likely model to account for the INS result, 
they proposed that (VO)2P207 has alternating antifer- 
romagnetic chains running in the b direction with unex- 
pectedly strong V-V exchange via the PO4 complexes. 
The exchange interaction Ji via the PO4 complexes is as 
large as 10~12 meV and the alternation ratio J2/J1 is 
estimated to be ~0.8. A possible importance of the su- 
perexchange pathways through the PO4 complexes was 
pointed out nearly a decade ago by Beltran-Porter et al. 
to account far the magnetic properties of various vanadyl 
phosphatesja and has recently been observed by INS in 
the dimer counterpart VODP04-iD20.Ll 

Another important but unresolved finding of Garrett 
et al. is the existence of a second excitation mode whose 
zone-center gap is 5.7 meV. They associated it with the 
two-magnon bound state as it lies just below 2A (A = 3.1 
meV being the gap for the one-magnon excitation deter- 
mined by Garrett et al). However, as they have pointed 
out, the substantial intensity of the second mode is diffi- 
cult to interpret for the two-magnon bound state. Frus- 
trated interchain coupling has recently been suggested to 
give rise to a triplet two-magnon bound state just below 
2A which is observable over a large part of the Brillouin 
zoncB However, the model is rather complicated and can- 
not reproduce correctly the intensity of the second mode 
especially near the zone center, where the distinct peak 
of the second mode is observed opposed to the theoretical 
prediction. 

In this paper, we report on the results of ^^P and 
Sly NMR on high-quality (VO)2P207 powder. We have 
found two independent magnetic subsystems existing in 
(VO)2P207 whose spin gaps are exactly the same as 
those of two distinct modes in the spin-excitation spec- 
trum. The existence of two magnetic subsystems with 
different spin gaps is supported by the recent high-field 
magnetization measurement. El From crystallographic con- 
siderations and a possible correlation between the ex- 
change interaction and the interatomic distance, we sug- 
gest that (VO)2P2 07 has two inequivalent alternating 
antiferromagnetic chains directed along the b axis. The 
role of the PO4 complexes as the exchange media and 
in determining the microscopic magnetic properties of 
(VO)2P207 will also be discussed in detail. 

This paper is organized as follows. We start with de- 
scribing experimental details in Sec. ^. The results of 
the NMR measurements will be presented in Sec. [H. 



Sec.0. 



The analysis of the line shape, the shifts and the nuclear 



spin- lattice relaxation rates will also be given in Sec. IH 



II. EXPERIMENTAL 



Polycrystalline sample of (VO)2P207 is prepared by 
following several steps of reaction. The first step is to 
synthesize VOP04-2H20 as a precursor by intercalating 
V2O5 with H3PO4 in ethanol at room temperature for 
about a month. The precursor is then heated in air at 
650 °C for 2 days to remove water completely. The final 
product (VO)2P207 is obtained by reducing oxygen from 
VOPO4 in Ar atmosphere at 700 °C for 2 days. The X- 
ray diffraction pattern of the obtained sample indicates 
that there is no trace of impurity phases. NMR measure- 
ments were performed with a standard phase coherent- 
type pulsed spectrometer. The ^^P NMR spectra were 
taken by integrating the spin-echo signal with a box-car 
averager while sweeping the external magnetic field at a 
fixed frequency. The ^^V NMR spectra were obtained by 
Fourier transforming the free-induction-decay signal af- 
ter a 7r/2 pulse at a fixed frequency and a magnetic field. 
Nuclear spin-lattice relaxation rates were measured by 
the saturation recovery method with a single saturation 
rf pulse for both the '^^P and ^^V sites. 



III. RESULTS AND ANALYSIS 

A. ^'^P NMR 



Section YV is devoted to the discussion of the experimen- 
tal results and their interpretations, followed by a con- 
clusion drawn from the present experimental findings in 



1. Line shape and the shifts 

The field-swept '^^P NMR spectra taken at various tem- 
peratures are shown in Fig. 0. Several peaks are found 
in the spectrum as expected from the crystal structure, 
and the line shape is strongly temperature dependent. 
At 4.2 K, a sharp single peak was observed near the zero 
shift K = 0, indicating the nonmagnetic ground state of 
(VO)2P207. The width of the spectrum is very narrow 
{^ 4 Oe) and is attributed to the broadening due to P nu- 
clear dipoles. The absence of inhomogeneous broadening 
of the spectrum assures the homogeneity of the sample 
in a microscopic scale and the extremely low concentra- 
tion of paramagnetic impurities. As the temperature in- 
creases, the spectrum splits into two peaks at first, and 
each peak further splits to show fine structure. Such a 
fine structure was not observed in the previous reports 
probably because of inhomogeneous broadening.ElEJ The 
multi-peak structure can be seen most clearly around 80 
K where the resonance lines are most distant from the 
K = position. There are four peaks visible in this 
temperature region, indicating that there exists at least 
four distinct P sites which feel different local magnetic 
fields from surrounding V atoms. Note that the num- 
ber of the observed peaks is less than the number of the 
crystallographically-inequivalent P sites. This suggests 



that the site distinction among some P sites is not clear 
and is indistinguishable in the present experiment. In 
the following, the observed peaks are referred to as the 
peaks 1 to 4 from low-field side to the higher. As the 
temperature is increased further, the fine structure is ob- 
scured and the spectrum merges into a single peak again 
(although the structure is visible at 250 K). 

The nonmonotonous variation of the line shape with 
temperature suggests the existence of several spin compo- 
nents having different spin susceptibility (xspin)- To trace 
the temperature (T) dependence of each spin component, 
we deconvoluted the spectrum into four gaussians and de- 
termined the shift of each peak Ki {i = 1, 2, 3, 4). The 
results are shown in Fig. g where Ki^s are plotted as a 
function of temperature. The T dependence of the Ki^s 
is distinct from each other in magnitude but is similar 
in that they all exhibit a broad maximum around|-8p K 
corresponding to the bulk magnetic susceptibility£3 At 
low temperatures, all the -ftT^'s decrease rapidly to almost 
zero owing to the gap in the spin-excitation spectrum. 

The difference of the shifts at the crystallographically- 
inequivalent nuclear sites is usually attributed to the dif- 
ference of the hyperfine coupling constants if there is only 
one kind of spin component. However, this is not the case 
in (VO)2P207. Figure ||(a) shows the shift of the peak 
3 {K3) plotted against the shift of the peak 2 {K2) with 
temperature the implicit parameter. If the T dependence 
of K2 and K^ is the same, the K3 versus K2 plot should 
lie on a straight line and its slope gives the ratio of the 
coupling constants A^/A2 {Ai being the hyperfine cou- 
pling constant at P sites associated with the peak i). It 
is obvious from Fig. ||(a) that K3 does not scale to K2 
at low temperatures (near the origin of the K3-K2 plot). 
There the K3 versus K2 plot concaves downward, which 
indicates that K2 increases more rapidly with increasing 
temperature than K3. This suggests that there are at 
least two independent spin components in (VO)2P207 
whose Xspin's are mutually different. 

It is interesting to note that at relatively high temper- 
atures, a nice linear relation is found between K3 and 
K2. This indicates that the T dependence of Xspin's re- 
sponsible for K3 and K2 is almost identical in this tem- 
perature region. This is not surprising because at tem- 
peratures much higher than the exchange interactions, 
spins behave paramagnetically and the asymptotic rela- 
tion Xspin ~ C/T with C the Curie constant for S = 1/2 
is expected to hold for all the spin components. From 
the slope of this linear portion, the ratio of the hyper- 
fine couphng constants A3/A2 is determined to be 0.88 
± 0.02. 

We have made a similar analysis for all possible pairs 
of Ki^s and have found that the observed four peaks can 
be divided into two groups. These are the peaks 1 and 
2 (group I) and the peaks 3 and 4 (group II) each of 
which shows distinct T dependence of the shifts. This 
indicates that there exist two independent spin compo- 
nents in (VO)2P207. Figure 0(b) shows the plots of Ki 
versus K2 and K^ versus K4^ similar to Fig. H(a). A 



linear relation is found between Ki {K3) and K2 {K4), 
indicating that their T dependences are identical aside 
from the T-independent prefactors (or Ai's). It is sur- 
prising that both plots lie on almost the same straight 
line in spite of the difference in the T dependence be- 
tween the groups. Indeed, the ratios of the coupling con- 
stants Ai/yl2 = 1.15±0.03and A3/A4 = 1.18±0.04 agree 
within experimental accuracies. This suggests that both 
P groups have similar local environment with regard to 
the site distinction among P sites within the same group. 
It is also noted that the intensity ratio of the two differ- 
ent P groups is close to unity (/1//11 — 1.1 ± 0.1). This 
indicates that the number of nuclei contributing to the 
NMR line of each P group is almost the same. 

As demonstrated in Fig. ^(a), the difference in the T 
dependence of the shifts at P sites belonging to the differ- 
ent groups is evident at low temperatures. This implies 
that the magnitude of the spin gap is different for the 
individual spin components. To estimate the gap of each 
spin component, we analyzed the Ki's based on the. sits-, 
ceptibility of a gapped one-dimensional spin system£3'Efl 
If the magnon dispersion along the chain is approximated 
by the quadratic form e{k) 2± A -|- c^fc^/2A near the bot- 
tom of the dispersion (fc = q — tt ~ 0), the T dependence 
of the susceptibility x in the low-temperature limit T <^ 
A is expressed as 
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where c is the spin velocity. We fitted the K2 (group 
I) and K3 (group II) data below 25 K to the form 
K = Xdia + aT-i/2exp(-A/T) where K^i^ is a T- 
indcpcndcnt diamagnetic contribution to the shifts. The 
obtained values are Xdia — —0.016 ± 0.005 % for both 
K2 and K3, a = 0.065 ± 0.002 K^/^ and A = 35 ± 2 K 
for K2, a = 0.079 ±0.002 K^/^ and A = 52±3 K for K3. 
The results of the fit are shown in Fig. || where we plot- 
ted {K - ATdia)!"^/^ against 1/T. Wc also plotted in the 
same figure the corresponding values of Ki and K4 with 
the appropriate scaling factors (i.e. the ratios A1/A2 and 
A3/A4 determined above) and -ftTjia. The correspondence 
between Ki {K3) and K2 {K4) is quite well. It should be 
noted that the estimated gap Aj = 35 K for the group 
I agrees within the experimental accuracy with the gap 
of the lower mode (36 K) observed in the INS on single 
crystals, while the gap An = 52 K for the group II is 
somewhat smaller than the gap for the upper mode (67 
K). 

The observed two independent spin components 
in (VO)2P207 may be associated with the two 
crystallographically-inequivalent magnetic chains run- 
ning along the b axis. The two chains have different 
alternation in the V-V distances, so that the gap may 
be different from each other. This would give a natural 
explanation for the existence of the second mode in the 
excitation spectrum, i.e., the second mode is not the two- 
magnon bound state but the independent magnetic exci- 
tation from another alternating antiferromagnetic chain 



with a larger spin gap. The fact that the NMR Unes 
for the groups I and II have almost the same intensity 
supports this interpretation. 

If we assume that the T dependence of Xspin's of the 
two spin components is identical at high temperatures 
which is justified by the presence of a linear region in 
the K^ versus K2 plot, the absolute values of the A^'s 
can be estimated from the slopes of the usual K-x plots. 
The plots of the KiS against x a-U have linear regions at 
high temperatures, the slopes of which yield Ax = 20 ± 2 
kOe/^B, ^2 = 17 ± 1 kOe//ZB, A3 = 15 ± 1 kOe/^e and 
A4 = 13 ± 1 kOe//iB- These values give the ratios A3/A2, 
A1/A2 and A3/A4 consistent with those determined in- 
dependently from the K-K scaling shown in Fig. ^. 



2. Nuclear spin-lattice relaxation rates 

The existence of two independent spin components 
with different energy gaps in (VO)2P207 will be demon- 
strated more directly if the nuclear spin-lattice relax- 
ation rates l/Ti for the different P groups can be mea- 
sured separately. However, the peaks overlap severely 
in the spectrum over a wide temperature range, which 
makes independent determination of 1/Ti rather difficult. 
We therefore measured 1/Ti in the limited temperature 
range from 12.5 to 30 K where the peaks for the differ- 
ent P groups are well articulated, and have successfully 
determined 1/Ti for the groups I and II. 

The recovery curves were taken at the peak position 
of each resonance line. To avoid excitations of nuclear 
spins belonging to the different group, we used a rela- 
tively small rf exciting field Hi. A typical strength of 
Hi was 14 Oe which exceeds the width of each peak and 
is still smaller than the line splitting (25-30 Oe at 31.4 
MHz where the measurements were performed). The re- 
covery curves can be fitted to an exponential function as 
illustrated in Fig. pi indicating that we reliably separate 
the signal from the different P group. 

The T dependences of 1/Ti are shown in Fig. M. It 
is clear that the 1/Ti for the different P groups show 
an activated behavior with different activation ener- 
gies. By fitting the data to an activation form l/Ti ex 
exp(— A'/T), we obtained the gaps Aj = 53 K and AJj 
= 71 K for the groups I and II, respectively. The gap 
is larger for the group II which is consistent with the 
shift measurement. It is also noted that the ratio of 
the two gap energies AJj/AJ = 1.34 ± 0.10 for 1/Ti 
agrees roughly with that of the shifts An/Aj — 1.48 
± 0.10. However, the absolute values of the gap deter- 
mined from l/Ti are larger than those obtained from the 
shifts asias-been observed in other gapped quantum-spin 
systems JIjLO 



B. ^W NMR 

The results of the '^^P NMR experiments show that 
(VO)2P207 has two independent spin components with 
different energy gaps for the spin excitations. This ac- 
counts for quite naturally the existence of two distinct 
excitation modes which are likely to come from the two 
inequivalent V chains running along the b axis. However, 
one of the gap energies estimated from the shift does not 
agree with the one observed in the INSu and magneti- 
zation measurementg3 that needs further examinations. 
This is partly because P nuclei feel transferred hyperfine 
fields from both spin components, which may result in 
the ambiguity in the estimation of the gap energy. The 
^^V NMR is free from such an ambiguity because of a 
dominant on-site hyperfine coupling and is expected to 
give a more precise measure of the T dependence of Xspin 
of each spin component. 

An intense free-induction-decay (FID) signal from V 
nuclei has been observed at low temperatures. The fre- 
quency spectrum obtained from the Fourier transform of 
the FID signal is isotropic with FWHM of about 70 kHz 
at 15.8 MHz and 4.2 K. The FWHM is nearly T inde- 
pendent and there is no indication of the line splitting up 
to 35 K. We could not trace the signal above 35 K due 
both to the poor signal-to-noise ratio and to the shorten- 
ing of the nuclear-spin relaxation time which prevents us 
from observing any transient signals by a pulsed NMR 
technique. 

Figure || shows the T dependence of the shift ^^K at 
the peak position of the spectrum. The shift is nearly 
constant below 10 K and varies strongly with tempera- 
ture above 10 K. In general, the shift K oi a. nucleus be- 
longing to the magnetic atom can be written as the sum 
of spin and orbital contributions, K — Ks-pm + ^orb- By 
assuming that Xspin vanishes at T = owing to the spin 
gap, the orbital shift K^rh at V sites is estimated to be 
0.212 ± 0.005 %. This is comparable with i^orb = 0.236 
± 0.005 % in the reference material V0HP0|iiiH20 with 
a similar local environment around V atoms .Ea The T de- 
pendence of i^spiu (= K — Ka^\^) is consistent with that of 
X given in Eq. (||), and the energy gap for the V sites is 
determined to be 68 ± 2 K by fitting Kg-pin to be propor- 
tional to X ill 6q. (y). This value is in excellent agreement 
with the gap of the second mode observed in the INS 
experiments. The observation of the gap of the second 
mode via NMR suggests that the second mode is not a 
branch for the two-magnon bound states but is simply a 
one-magnon propagating mode. 

The inset of Fig. [s] is a plot of ^^if versus bulk x with 
temperature the implicit parameter. It is obvious from 
the figure that the ^^K does not scale to the bulk x- The 
^^K-x plot concaves upward, which indicates that Xspin 
responsible for the observed ^^K decreases more rapidly 
with decreasing temperature than the bulk x- The result 
of the ^^K-x plot also suggests that the bulk x includes 
an additional contribution from the spin component with 



a gap smaller than that observed via ^"'^V NMR. Unfor- 
tunately, we could not detect a transient V signal cor- 
responding to the spin component with a smaller gap 
energy, which results probably from the short nuclear 
spin-spin relaxation time T2 of such V sites. 

The T dependence of 1/Ti measured for the observed 
FID signal is shown in Fig. ||. An activated behavior with 
a gap of 66 ± 2 K has been observed down to 14 K, be- 
low which the 1/Ti starts to deviate from the activation 
law probably because of extra relaxation processes. The 
obtained gap for 1/Ti agrees quite well with that deter- 
mined from thC|-^^-ftr measurement, which is expected in 
the low-T limit .Ej This is contrasted with the result at P 
sites where the gap energies for the shifts and l/Ti's are 
different. The reason for this difference between V and 
P sites is not clear, although the different wave-vector 
dependence of the hyperfine form factor is likely to be 
the case. 



IV. DISCUSSION 

It is now evident from our NMR experiments that 
(VO)2P207 is composed of two independent magnetic 
subsystems with different spin-gap energies. Since 
the different magnetic properties of the system gener- 
ally come from the different crystallographic environ- 
ment, it is reasonable to associate them with the two 
crystallographically-inequivalent V chains running along 
the b axis which (VO)2P207 inherently possesses. The 
two-chain model gives consistent explanations for most of 
the present experimental results such as the observation 
of two distinct spin gaps (and Xspin's), the intensity ratio 
of the resonance lines for the two different P groups, and 
the underestimation of one of the spin-gap energies. This 
model also accounts for the existence of the two distinct 
magnetic modes, which have been thought to be a one- 
magnon and the two-magnon bound states of the same 
alternating antiferromagnetic chain. 



A. Qualitative aspects 

Before discussing the experimental results in detail, let 
us review the crystal structure of (VO)2P207. As has 
been noted previously, (VO)2P207 has eight inequiva- 
lent V and P sites, and 36 oxygen sites in the unit cell. 
This makes a crystallographic difference between the two 
kinds of V chains running along the b axis. The chain 
"A" consists of V atoms on V1-V4 sites andr-the chain 
"B" includes V5-V8 sites as shown in Fig. |i|.EI The V- 
V distances dyy's in both the chains A and B alternate 
to give rise to the alternations of the exchange interac- 
tions, but the alternations in dyv's are slightly different 
between the two chains. 

The difference of the spin gap and the resultant differ- 
ence of Xspin between the chains can be understood qual- 



itatively by considering that the different alternations in 
rfvv results in the different alternations in the exchange 
interaction as well. This is very likely because the ex- 
change interaction is expected to depend sensitively on 
the interatomic distance through the overlap of the re- 
lated electronic orbitals. In the alternating antiferromag- 
netic chains, the shape of the X"^^ curve varies strongly 
with the ratio a = J2/ JiXl^ Therefore, if the chains A 
and B have different values of a owing to the different 
alternations in dyv; they exhibit different T dependence 
of Xspin's. The difference of the spin gap can also be 
explained because the spin gap of the alternating chain 
depends on the ratio a as well as the magnitude of Ji 
and J2- 

We did not mention explicitly the exchange pathways 
in the above argument. Garrett et al. have pointed out 
that the exchange via the PO4 complexes in the b direc- 
tion is essential in determining the magnetic properties 
of (VO)2P207 as alternating antiferromagnetic chains.cl 
Most of the present experimental results are consistent 
with such exchange pathways and are compatible with 
the two-chain model as described below. 

For the magnetic connection along the b axis, two PO4 
complexes in between the VO5 pyramids should be the 
important exchange pathways. In the chain A, the ex- 
change between VI and V4 (V2 and V3) sites is mediated 
by the PO4 complexes containing P6 and P8 (P5 and P7) 
sites. It includes the distant V-0-P-O-V paths. There 
is also the exchange coupling between V4 and V2 (V3 
and VI) sites via the shared VO5 pyramidal edge hav- 
ing shorter V-O-V paths. This is illustrated in Fig. nO 
together with the exchange paths for the chain B. As a 
natural consequence of the exchange via PO4 along the 
b axis, eight inequivalent P sites are divided into two 
groups depending on to which chain they belong as the 
exchange paths: P5-P8 sites belonging to the chain A 
and P1-P4 sites belonging to the chain B. This is consis- 
tent with the observation of two different P groups via 
■^^P NMR. In addition, since the number of the sites is 
the same for both P groups in this model, the intensity 
of the NMR lines should be the same. This is what was 
observed in the '^^P NMR experiments. 

The exchange coupling via the PO4 complexes explains 
why the susceptibility of each spin component can be 
traced selectively by means of the '^-'^P NMR. In the two- 
chain model, the shift Kp^ at P sites in the chain A can 
be written as A a = flAAxfpin + aABxfpin- Here x^i„ and 
xHjin the spin susceptibilities of the chains A and B; aaa 
and OAB the coupling constants at P sites in the chain A 
for the hyperfine fields transferred from the V sites in the 
chains A and B, respectively. Similarly, for P sites in the 
chain B, the shift K^ can be written as A'b = '^BAX^in"'" 
flBBX^in with Aba and abb defined in a manner similar 
to AAA and Cab- Since the transferred hyperfine field at 
nuclear sites in between the exchange-coupled magnetic 
atoms arises from the electron-spin transfer similar to 
that for the exchange, the hyperfine field at the P site 



is expected to come mainly from the V atoms in the 
same chain. This means that IaabI and |aBA| are much 
smaller than JcaaI and IobbI, indicating K{^ and K^ to 
be roughly proportional to x^in ^'^'^ ^?pin' respectively. 

It is expected that the interchain exchange coupling 
in the c direction is weakened by the presence of the 
intrachain coupling via the PO4 complexes. This is be- 
cause molecular orbitals of the PO4 complex should be 
configured to connect the neighboring V atoms in the h 
direction, which may be a microscopic origin of the domi- 
nant hypcrfinc coupling between the P nucleus and the V 
moments in the same chain. As expected, the magnetic- 
excitation spectrum of (VO)2S'207 exhibits a very weak 
dispersion in the c direction.Q We notice here that the 
weak ferromagnetic coupling in the a direction acts be- 
tween the same kind of chain so that the chains A and B 
are magnetically decoupled. It is to be emphasized that 
the periodicity of the lattice is unchanged by the resul- 
tant grouping of V and P atoms: the magnetic unit cell is 
identical with the chemical one for both magnetic chains. 
Therefore, the spin- wave dispersion is observed in exactly 
the same Brillouin zone. This is consistent with the INS 
results. 

Although the contribution from the different kind of 
chain to the hyperfine fields at P sites is considered to be 
small, it affects the estimation of the gap energies from 
Kp^ and K-q. This is serious for the P group belonging 
to the chain with a larger spin gap for the following rea- 
son. Suppose that the spin gap Ab of the chain B is 
smaller than the spin gap Aa of the chain A. Then xfpin 
starts to grow at lower temperatures than x^in ^^'^ there 
may be a temperature region in which xSjin grows rapidly 
with increasing temperatures while x^in i^ ^^ill vanish- 
ingly small. In such a temperature region, K^ may be 
approximated by aBBxfpin because both oab and x^in 
are small. Hence the value of Ab will be estimated cor- 
rectly from K-Q- On the other hand, for P sites in the 
chain A, the contribution from the different kind of chain 
flABxfpin rnay be comparable with, or greater than that 
from the same chain OAAX^in- ^^ ^ result, Kp^ would 
start to increase at lower temperatures than x^in- This 
is inevitable as far as oab is finite andpAA > Ab, and 
will result in the underestimation of Aa-E3 In the present 
NMR experiment at P sites, the gap An of the group II 
is smaller than the gap of the upper mode in the exci- 
tation spectrum, while Aj agrees well with the gap of 
the lower mode. In the light of the above argument, the 
underestimation of An is considered to result from the 
small but finite contribution from V sites in the different 
kind of chain with a smaller gap. At this point, we may 
associate the P sites belonging to the group II (I) with 
the chain having a larger (smaller) spin gap. The value 
of the gap determined from the ^^V NMR experiment 
agrees well with the INS result and is considered to be 
more accurate for the larger spin-gap energy. 



B. Quantitative aspects 

Most of the present experimental results can be ex- 
plained qualitatively in terms of the two-chain model by 
taking account of the relatively strong V-V exchange via 
PO4 in the b direction. The model is consistent also with 
the INS result and therefore seems to be appropriate for 
the description of the magnetic behavior of (VO)2P207. 
However, one crucial question arises here whether the 
small differences in dvv can make the spin gaps of the 
two alternating chains be different by a factor of nearly 
two. To answer this question, we estimated the exchange 
interactions in each alternating chain. It is found that the 
difference of the exchange is not very large and may be 
accounted for by the strong dependenc&jof-the exchange 
interaction on the interatomic distance£jc3 

For the analysis, we take the values of the spin gap to 
be 35 K and 68 K. The former is the one determined from 
the shift at P sites belonging to the group I, and the latter 
from the shift at V sites. Both are in excellent agreement 
with the zone-center gaps for the lower and upper modes 
in the spin-excitation spectrum, respectively. The zone- 
boundary energy ezB taken from the spin-wave dispersion 
along the b axis is 15.4 meV (a? 179 K) and is assumed 
to be the same for both chains.Q Then J = ( Ji -I- J2)/2 is 
calculated from the relation tt J/2 = ezB to give the ratio 
A/ J, and the alternation parameter 5 = (1 — a)/(l-|-a) is 
determined by using tiue 5 dependence of A/ J presented 
by Uhrig and Schluz.EZi From the values of J and ^, the 
exchange interactions Ji and J2 are calculated. The re- 
sults are summarized in Table ||. 

We have assigned the chain A to have a larger exchange 
alternation (i.e., a larger value of 5) resulting in a larger 
spin gap. This is based on the following two assump- 
tions: 1) the exchange via the PO4 complexes is stronger 
than that through the shared VO5 pyramidal edge, 2) 
the interaction is weaker for longer dyv for both the ex- 
change paths. The V-V distance via the PO4 complexes 
in the chain A is ~5.14_4, which is shorter than that in 
the chain B of ^5.16 A.c3 The largest exchange coupling 
(Ji = 136 K) is therefore assigned to the V-0-P-O-V 
paths in the chain A. The V-O-V path in the chain A 
is then assigned automatically to have the smallest cou- 
pling constant (J2 = 92 K) without any additional as- 
sumptions. The V-V distance is ~0.02 A longer than the 
corresponding dyv in the chain B of ~3.21 A and hence 
the assignment is self-consistent. Even if we regard the 
shared pyramidal edge as a dominant exchange pathway, 
the assignment can be made in the same way without in- 
consistency and now the spin gap of the chain A should 
be smaller than that of the chain B. 

It appeared that the exchange interactions of the chain 
A should differ from those of the chain B by ^10 K ('-^lO 
% in relative strengths) to give rise to the observed dif- 
ference of the spin gap between the two chains. Since the 
two-chain model is based on the crystallographic inequiv- 
alence of the chains, the difference of the exchange should 



be ascribed to the small differences (~0.02 A for both 
the exchange paths) in dyv ■ With regard to the V-V ex- 
change via the shared pyramidal edge, the difference may 
be attainable because in some copper salts with edge- 
sharing Cu-0 bondings, the 0.02 A difference in the Cu- 
Cu distance-gives rise to the difference in the exchange 
by ~100 K.E3 The strong dependence of the exchange, 
the interatomic distance is also known empirically,Ej 
putting the basis on the dependence of overlap integrals 
on distance. For V-V pairs sharing common pyramidal 
edge, the empirical law predicts the-ftxchange to decrease 
with increasing dyv as {dyv)~^^E3 The application to 
(VO)2P207 yields ^^6 % decrease of the V-V exchange 
on increasing dyv from 3.21 to 3.23 A, which is com- 
parable with the above estimation. Although such an 
empirical relation is not known for V-V pairs connected 
by the PO4 complexes because of the complicated and 
unresolved interaction pathways, it is plausible that the 
exchange is also sensitive to dyv as the exchange via the 
common VO5 edge. The small difference in dyy's be- 
tween the two crystallographically-inequivalent chains is 
therefore considered to be crucial in making them mag- 
netically distinct from each other. Quantum chemical 
calculation of the exchange interaction via the PO4 com- 
plexes is desirable for more thorough understanding of 
the magnetic properties of (VO)2P207. 



V. CONCLUSION 

We have presented the results of the extensive NMR 
study on the quasi one-dimensional gapped quantum an- 
tiferromagnet (VO)2P207. As revealed by the detailed 
line-shape analysis of the ^^P NMR spectra and the mea- 
surement of 1/Ti at ^^P sites, (VO)2P207 has two inde- 
pendent spin components with different gap energies for 
the spin excitations. The gap energies are determined 
from the temperature dependence of the shifts at '^^P 
and ^^V sites to be 35 K and 68 K, and are found to 
coincide with those observed in the INS experiments. 

The agreement of the gap energies between NMR and 
INS experiments suggests that the two magnetic modes in 
the excitation spectrum arise from the two distinct mag- 
netic subsystems which (VO)2P207 inherently possesses. 
We proposed that two crystallographically-inequivalent 
V chains running along the b axis are magnetically dis- 
tinct to contribute independently to the spin-excitation 
spectrum. The two chains have different exchange alter- 
nations and spin gaps, which may result from the small 
difference in the V-V distances between the chains. It 
is also suggested that the PO4 complexes act as the im- 
portant superexchange pathways. This is essential in di- 
viding the P sites into two groups and in observing the 
susceptibility of each alternating chain selectively via '^^P 
NMR. Quantitative estimates of the exchange interac- 
tions show that the difference of the exchange is not very 
large and may be attributed to the difference in the V-V 



distances between the two chains. 

To conclude, (VO)2P207 is a complicated magnetic 
system containing two kinds of antiferromagnetic linear 
chains with different alternations in the exchange inter- 
action. The energy for the spin excitations is different 
in each alternating chain, which is the origin of the two 
distinct magnetic modes in the excitation spectrum. 
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FIG. 1. Polyhedral view of the crystal structure of 

(VO)2P207. Also indicated are site indices of the inequivalent 
V and P atoms included in each polyhedron. The structural 
units indicated are stacked in the a direction. See text for the 
labeling of the chains directed along the h axis. 



FIG. 2. Temperature variation of the P NMR spectrum 
taken at 22.6 MHz. Dashed line indicates the position of the 
zero shift {K = 0) for "'^P nuclei calculated using the nuclear 
gyromagnetic ratio '*^7n ~ 17.235 MHz/T. Numbering of the 
peaks (see text) is shown for the spectrum at 77.4 K. 
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FIG. 3. Temperature dependence of the shifts at "^^P sites. 
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FIG. 5. Plots of {K - Ka\^)T^^'^ against l/T. The data 
for K\ and A'4 are plotted by multiplying the scaling factors 
A2/A1 and A4/A3 to absorb the difference of the hyperfine 
coupling constants within the groups. The lines indicate the 
activation law {K - A'dia)T^''^ oc exp(-A/T) with A of 35 K 
and 52 K for the groups I and II, respectively. 
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FIG. 4. (a) The shift of the peak 3 {Kz) plotted against 
the shift of the peak 2 {K2) with temperature the implicit 
parameter. The straight line is a fit of the data from 100 to 
250 K. (b) Plots of Ki versus K2 and K3 versus Ka. The 
line is a guide to the eyes of which slope corresponding to the 
average of A1/A2 and Az/A^. 
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FIG. 6. Typical examples of the recovery curve of "^^P 
nuclear magnetization. The data were taken at 31.4 MHz 
and 20 K. Inset shows the field-swept ^^P NMR spectrum 
at the same frequency and temperature. Arrows in the inset 
indicate the positions of the external fields where the recovery 
curves are measured. 
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FIG. 7. Temperature dependences of l/Ti at the ^^P 
sites belonging to the groups I and II. Solid lines are the 
results of the least-squre fit of the data to the activation law 
l/Ti oc exp(- A'/T) with A' of 53 K and 71 K for the groups 
I and II, respectively. 
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FIG. 8. Temperature dependence of the shift ^^K at the 
^^V sites. Solid line is a fit of the data with a gap of 68 K 
(see text). Inset shows ^^ A' plotted against the bulk magnetic 
susceptibility. 
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FIG. 10. Schematic description of the exchange interac- 
tions in (VO)2P207. The assignment of the inequivalent P 
sites to each alternating chain is also shown. 



TABLE I. 

(VO)2P207 


Spin gaps and exchange interactions in 


A (K) Ji (K) J2 (K) J2/J1 


chain A 
chain B 


68 136 92 0.67 
35 124 103 0.83 
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